Propagation of Nonclassical Radiation through a Semiconductor Slab 
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Based on a microscopic derivation of the emission spectra of a bulk semiconductor we arrive at a 
clear physical interpretation of the noise current operators in macroscopic quantum electrodynamics. 
This opens the possibility to study medium effects on nonclassical radiation propagating through 
an absorbing or amplifying semiconductor. As an example, the propagation of an incident squeezed 
vacuum is analyzed. 

PACS numbers: 42.50.Nn, 78.20.-e, 42.50.Ct 



The progress of experimental techniques has rendered 
it possible to almost completely measure the quantum 
statistics of radiation-matter systems, so that correla- 
tion functions up to, in principle, arbitrarily high orders 
can be detected (for a review of the various methods, 
see, e.g., Ref. [H). Within the framework of quantum 
optics, the theoretical study of the interaction of light 
with complex material systems is typically based on sim- 
plifying model systems, cffectivc-Hamiltonian schemes or 
related semi-phenomcnological concepts rather than rig- 
orous microscopic calculations 0. 

On the other hand, in many-particle quantum theory, 
including semiconductor theory, much effort has been 
spent on the development of methods for microscopically 
describing complex material systems. This includes the 
description of coherent optical interactions by semicon- 
ductor Bloch equations and semiconductor luminescence 
equations as well as the development of noncquilibrium 
Green's function methods 0, Q ■ In principle, these meth- 
ods lead to infinite hierarchies of correlations, which are 
usually treated by properly developed methods of trun- 
cations and/or decorrelations. 

It is well known that semiconductors can be used to 
generate nonclassical radiation [Bj]. In particular, the de- 
velopment of nano-structured systems has opened new 
possibilities of the generation and application of non- 
classical radiation in integrated systems. For example, 
the correlated emission of singlephotons can be demon- 
strated by using quantum dots [6(] and bound excitons in 
semiconductors Experiments with quantum wells [§] 
and quantum dots 0, [l(| EH also show the potential of 
semiconductors for the generation of entangled photons, 
which are of interest in quantum information processing. 

In order to properly describe the generation and/or 
propagation of nonclassical radiation through complex 
material systems such as semiconductor slabs, one has 
simultaneously to deal with both higher-order radiation- 
field correlation functions and many-particle quantum 
statistics of the material system. Within the frame of 
macroscopic quantum electrodynamics (QED), methods 



of describing the quantized electromagnetic field in lin- 
early responding (dispersing and absorbing/amplifying) 
media have been developed, wit h sp ecial emphasis on 
quantum- noise effects 12, 13, 14, TH, 3, 13] ■ Based on 



given constitutive relations of the material system, field 
correlation functions of arbitrary order can be calculated 
in this way. On the other hand, methods of many-particle 
theory can be used to treat the radiation-matter interac- 
tion within the frame of microscopic QED, thereby cal- 
culating the relevant electromagnetic properties of the 
material system. Such methods have been used to study 
the propagation of radiation in semiconducting material, 
including amplification and lasing T^, 2^], with the re- 
striction to the intensity and the emission spectrum of 
the radiation, because of the elaborateness of the prob- 
lem. 

In the present paper we make the attempt to close the 
gap between macroscopic and microscopic QED, lead- 
ing to a clear physical interpretation of the noise-current 
density introduced in macroscopic QED. In particular, 
appropriate combination of the results of the two meth- 
ods offers the possibility to calculate arbitrary radiation- 
field correlation functions, as we demonstrate by study- 
ing the propagation of squeezed light through an absorb- 
ing/amplifying semiconductor slab. 

We begin with the operator of the vector potential in 
Coulomb gauge, A(r, t), which obeys the inhomogeneous 
wave equation 



□ A(r,i) = -Moj(r,t) 
and the equal-time commutation rule 



A(r,i)^A(r',t) 



ihc 



5 T (r-r'), 



(1) 



(2) 



where <5x(r) is the transverse tensorial delta function. 
The operator of the transverse current density, j(r,i), 
can be decomposed into the current density j m cd(r, t) as- 
sociated with the medium under consideration and an 
auxiliary, externally controlled c-number current density 



2 



jext(r,t) which is commonly set zero at the end of the 
calculations. 

Starting from microscopic QED and employing Green's 
function (GF) techniques, for a dielectric medium in the 
steady state, the expectation value of the vector potential 
(in the Fourier domain) reads as 

(i(r, uj))= - M o / dV V rot (r, r', u>) • j cxt (r', u), (3) 



where, in linear approximation, the inverse of the re- 
tarded photon GF D rct (r, r', ui) is given by 



D> 



l (r,r',uj) = 05(r-r')- P rct (r, r', w), (4) 



with the retarded polarization tensor P rot (r, r',u>) being 
related to the complex (dielectric) susceptibility tensor 
X (r, r', uj) — x '( r j r 'i w ) + i X "( r > r '; w ) 01 the medium 
according to 



P ret (r, t',lj) 



X*(r,r',u;) 



(5) 



The Keldysh components of the photon GF, 
*D>(r,r',t-t') = V<(r',r,t' -t) 

r (A(r,i)A(r',t'))-{A(r,t))(A(r',0)j, (6) 

describe field-field fluctuations. Solving the Dyson equa- 
tion for these components, one obtains the optical theo- 
rem 



J d 3 n J d 3 r 2 D^{T, ri ,u) 



x P<(r 1 ,r 2)W )D I f'(r a ,r>) ) (7) 



where P^ are the corresponding Keldysh components of 
the polarization tensor, and summation over repeated in- 
dices is understood. Within the framework of macro- 
scopic QED, Eq. JT]) is regarded as being the operator- 
valued inhomogeneous wave equation that corresponds 
to the Maxwell equations of the transverse part of the 
macroscopic electromagnetic field in a linear dielectric 
medium. Hence, j m ed(r,w) must have the form 

j mcd (r,u)=-e uj 2 / dV^(r,r',w)A(r',w) +j N (r,w), 



(8) 

and in place of Eq. ([3]) we obtain the operator-valued 
equation 

A(r,w) 

=-Mo / d 3 r' < D rct (r,r',^).[j N (r',^)+j cxt (r', W )l , (9) 



where the noise current density operator Jn(t, uj) obeys 
the commutation relation 11611 



[jN(r, — lo) = jjj(r, uj)], which ensures the validity of 
Eq. ©. 

With the analog of the optical theorem ([7]) in macro- 
scopic QED, the relation 



P>{v,v',t-t') = P<{v',v,t'-t) 

= g(jN(r,t)j N (r',t')) 



(11) 



can be derived. Moreover, using Green's function tech- 
niques, a Bethc-Salpeter equation between the polariza- 
tion function and the correlation function of the medium 
current density can be derived [l8| . In this way, the fluc- 
tuation of the noise current density in macroscopic QED 
can be directly related to the fluctuation of the micro- 
scopically well-defined and observable medium current 
density. 

In the following we will deal with the propagation of 
TE-polarized radiation along the x axis, through a semi- 
conductor slab of thickness L which is infinitely extended 
in the y-z-planc. For simplicity, we assume that the elec- 
tric field is polarized along the y-axis A=(0, A, 0). Then, 
neglecting spatial dispersion, the complex refractive in- 
dex inside the medium, n = n'+in" , is obtained from 
n 2 H = l + x(w). 

Using Green's function technique, the spontaneous 
emission of the slab is studies in Ref. [3, A gen- 
eralization of the results including spatial dispersion and 
providing exact relations between (spontaneous and stim- 
ulated) emission and linear absorption can be found in 
Ref. [22j. The intensity of the (spontaneously) emitted 
radiation can be given by 



= &(w)P(w)A>(u;) , 



(12) 



where Dq is the vacuum-induced contribution to the 
photon spectral function, P = 2iImP rct is, accord- 
ing to ((5]), related to x" associated with linear absorp- 
tion/amplification, and b is defined as the ratio between 
the recombination rate P < of electron-hole pairs and P. 
The function b{uS) characterizes globally (i.e., inside and 
outside the slab) the emitted radiation and, as such, b is 
accessible to direct observation in experiments. It gen- 
eralizes Planck's formula for the black body radiation to 
nonequilibrium radiation of an excited medium in the 
steady state. 

On the basis of Eq. ([9]), input-output relations can 
be derived, by introducing bosonic quasiparticle annihi- 
lation operators (normalization A/±) 



fL/2 

c±(uj)=j\f±(uj) dx [t 

J-L/2 ^ 



ox'/c ± e -™"*'/Aj N ( x > rU) ) 

(13) 

expressed in terms of the noise current density op- 
erator and associated with the slab-radiation excita- 
13ll2lj]. On the other hand, the poles of L> rct yield 



tions 



[j]sr(r, ui), j Jj (r', uj')] — 2eohuj 2 '^x > "(r, r'u))S(u— uj') (10) the polaritonic dispersion relations for the slab, which 
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allows one to relate the operators c± to the polaritonic 
annihilation operators [23| . Assuming that the input field 
is in the vacuum state, the intensity of the output field 
is obtained in just the same form as in Eq. (fT2")> . by iden- 
tifying (c 1 (u>)c(u)) = b(uj), where c(oj) = c + (uj) + c_(lj). 

Steadily excited semiconductors in quasi-equilibrium 
are of particular interest, such as exciton gases gener- 
ated at low up to moderate excitation and light-emitting 
diodes working at high excitation. For quasi-equilibrium, 
due to the Kubo-Martin-Schwinger condition (3], the 
function b(uo) develops into a Bose distribution b(uj) = 
(exp [f3(huj — — . The chemical potential /i starts 
at fx = for complete thermal equilibrium and charac- 
terizes the degree of excitation beyond the thermal one 
for /i > 0. The crossover from absorption (x" > 0) to 
gain < 0) appears at huj = fi. By expanding the 
product 6(w)x"(w) in Eq. (fT2|) at Huj = fi, it is seen that 
the spontaneous emission remains finite at the crossover, 
it is given by the slope of the absorption function x" ■ 
Since both x" an d b switch their signs, the spontaneous 
emission stays positive in the whole frequency region, as 
it should be. Thus we obtain a unified description of 
absorption and amplification in semiconductors. 

In order to demonstrate the usefulness of the estab- 
lished interrelation of microscopic and macroscopic QED, 
we consider the propagation of squeezed light through a 
semiconductor slab. Consider the squeezed vacuum state 



T=3K 
T=300K 



|V>)sv = S\lp) v , 



where the squeeze operator, 



(14) 



S = exp I J dio [£* a(oJo+oj)a(uJo— u>) — h.c] j (15) 

acts on the multimode vacuum state The squeeze 

parameter is £ = |£|e l< ^ 5 , and a(ui) are the photonic an- 
nihilation operators of the incident field on the left-hand 
side of the semiconductor slab. The input field on the 
right-hand side is assumed to be in the vacuum state. 

Using Eqs ©, (JH]), and (Hi]), the normally ordered 
squeezing spectrum 5 sq = J cLuj'u!q(: A(x, uj)A(x' , ui') :) sv 
is derived for x~x'~L as 

^ q = /H + -^L{|TM| 2 sinh 2 |e| 

+ e-^L/c^, [T»] 2 cosh|£| sinh|e|+c.c.}. (16) 

Here we have assumed that the detector is placed on 
the right-hand side of the slab. The total transmission 
coefficient of the slab, T(lo), reads in terms of the internal 
reflectivity r{uj)=[l-n(uj)]/[l+n{u)]e lu " l ^ L / c as [25| 



4n(w) 



Juj(n(ui)-l)L/c 



(17) 



[l+n(iu)} 2 [l-r 2 (uj)} 

From Eqs IT2|) and (fT6|) it is evident that the behavior of 
the function S sq strongly depends on the behavior of the 
denominator 1 — r 2 (oj). 




hio — E g meV 



FIG. 1: Output squeezing spectrum for a GaAs slab (L — 
25fim) near the ls-excitonic resonance (indicated by an ar- 
row) of bandwidth = 0.2 meV. The squeezing spectrum is 
shown for two temperatures and for |£| = 0.2 (a), 1.2 (b). 
The dashed lines indicate the maximum and minimum noise 
level of the squeezed input field. 



So, the squeezing spectrum shows an oscillating behav- 
ior with the maxima at the frequencies ui s = L ^J^ s with 
integer s (Fabry-Perot resonances). 

In Fig. [T] we show the influence of a single excitonic res- 
onance at E x =fijjjQ on squeezed light propagating through 
a semiconductor slab, for different values of the squeezing 
strength |£| where the Lorentz oscillator model has been 
used to simulate excitonic absorption [26| . The input 
(squeezed white noise) spectrum is depicted with dashed 
lines. The upper (gray) and lower (black) curves repre- 
sent the maximum and the minimum noise level of the 
squeezing spectrum of the radiation after transmission 
through the slab. For high temperatures, the emission 
spectrum I (to) plays a more pronounced role and it leads 
to a significant decrease of the nonclassical properties of 
the input field. 

The squeezing spectrum for a highly excited semicon- 
ductor is shown in Fig. [2^a) for |£| = 0.2. The suscepti- 
bility function has been calculated in the effective pair- 
equation approximation [26| and is given in Fig. [2Kb) . In 
the output field, the squeezing is nearly destroyed, due 
to the dominance of the incoherent emission in the gain 



4 





30 


CO 

+^ 

'3 


OA 














cr 


10 












(a) 




4 6 8 10 

ftuj — E„ meV 



a3 



.(b) 




p X " ~ 


\. 2 


4 


J 6 8 10 






/ hu — E g meV 









FIG. 2: Squeezing spectrum after transmission of a GaAs 
slab at temperature T=3K, pumped to a carrier density of 
n = 3 • 10 16 cm~ 3 (a) and the corresponding imaginary part of 
susceptibility (b). The input (dashed) and transmitted (solid) 
squeezing spectra are shown for |£| = 0.2. An arrow indicates 
the chemical potential fi and the inset shows its vicinity. 



region. Minor squeezing is only preserved in the vicinity 
of the crossover at Huj = fx, where x" = 0. 

In conclusion, we have studied the relation between mi- 
croscopic and macroscopic QED, leading to a clear physi- 
cal interpretation of some basic quantities. By combining 
the input-output formalism of macroscopic QED with 
results of microscopic QED, the treatment of the prop- 
agation of nonclassical light in complex media becomes 
possible. As an example, we have studied squeezed-light 
propagation through a semiconductor slab. 

This work was supported by the Deutsche Forschungs- 
gemeinschaft, Sonderforschungsbereich 652. 
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